Short-Time Operator Product Expansion 
for rf Spectroscopy of a Strongly-interacting Fermi Gas 
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Universal relations that hold for any state provide powerful constraints on systems consisting of 
fermions with two spin states interacting with a large scattering length. In radio-frequency (rf) 
spectroscopy, the mean shift in the rf frequency and the large frequency tail of the rf transition rate 
are proportional to the contact, which measures the density of pairs with small separations. We 
show that these universal relations can be derived and extended by using the short-time operator 
product expansion of quantum field theory. This is a general method for identifying aspects of 
many-body physics that are controlled by few-body physics. 
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Trapped ultracold atoms allow the study of the few- 
body physics and the many-body physics of systems in 
which the fundamental interactions between the con- 
stituents are understood and can be controlled exper- 
imentally. The possibility of making the interactions 
between the atoms arbitrarily strong presents a chal- 
lenge to theory, because many theoretical methods break 
down when interactions become too strong. Specifying 
the interactions between the constituents of the many- 
body system in terms of 2-body interactions is a trivial 
example of a connection between few-body physics and 
many-body physics, but there can be deeper connections. 
Among the simplest examples of strongly-interacting sys- 
tems are ones that consist of atoms that interact with 
a large scattering length [lj. The 2-body problem for 
such atoms can be solved analytically, the 3-body prob- 
lem can be solved exactly numerically, and the 4-body 
problem is becoming tractable using modern computers. 
This makes it an interesting system for studying how non- 
trivial aspects of many-body physics can be controlled by 
few-body physics. 

A powerful experimental tool for studying ultracold 
atoms is radio frequency (rf) spectroscopy, in which 
atoms in one hyperfine spin state are excited into a dif- 
ferent spin state. Pioneering applications of this method 
were measurements of the binding energy of weakly- 
bound diatomic molecules of 40 K atoms @ and the study 
of the pairing gap in a many-body system of 6 Li atoms 
The method had been extended to allow the spa- 
cial resolution of trapped systems [ij and the momentum 
resolution of the excited atoms [f| • A recent review of rf 
spectroscopy has explored the relation to photoemission 
experiments on high-temperature superconductors @. 

Universal relations between various properties of a 
system that must be satisfied in any state can provide 
powerful constraints on theoretical methods. Shina Tan 
has derived universal relations for systems consisting of 
fermions with two spin states that interact with a large 



scattering length a Q ■ These relations all involve a prop- 
erty of the system called the contact. It is an extensive 
quantity that can be expressed as the integral over space 
of the contact density, which is proportional to the num- 
ber of pairs with different spins per (volume) 4 / 3 @|. The 
Tan relations include the coefficient of the 1/k 4 tail in 
the momentum distributions at large momentum, a de- 
composition of the total energy into terms that are in- 
sensitive to short distances, the rate of change of the free 
energy due to a change in a, the relation between the 
pressure and the energy density in a homogeneous sys- 
tem, and the virial theorem for a system in a harmonic 
trapping potential Tan derived his relations within 
the framework of the many-body Schrddinger equation 
using novel methods involving generalized functions. In 
Ref. [8| , Braaten and Platter rederived the Tan relations 
within the framework of quantum field theory using stan- 
dard renormalization methods together with the short- 
distance operator product expansion. The Tan relations 
have also been rederived using less formal methods [9(. 

The contact also plays an important role in the radio- 
frequency (rf) spectroscopy of ultracold atoms. In a 
many-body system consisting of atoms in spin states 1 
and 2, an rf signal with frequency oj can be used to trans- 
fer atoms from the spin state 2 into a third spin state 3. 
There are sum rules that constrain the rf transition rate 
r(w) Of particular interest is the case of fermionic 
atoms for which every pair has strong interactions de- 
termined by large scattering lengths eti2, ai3, and 023. 
An example is the lowest three hyperfine spin states of 
6 Li atoms at generic values of the magnetic field. In this 
case, the sum rules are 



duj r(cj) = Trn 2 N 2 , 



(la) 



dujujT(uj) = (n 2 /Am)(a^ ~ a^)C 12 , (lb) 



where Q is the Rabi frequency determined by the 
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strength of the rf signal, N 2 is the number of atoms of 
type 2, and C12 is the contact for atoms of types 1 and 2. 
We measure the rf frequency u relative to the hypcrfinc 
frequency difference between atoms 2 and 3, and we set 

h 



and antiquark at a point with total energy E and with 
equal and opposite momenta. They are subsequently 
transformed by the strong interactions of QCD into one 
or more hadrons. The inclusive cross section u(E) for 



1. The sum rule in Eq. (jlbl 
mean shift in the rf frequency due to interactions is pro- 
portional to the contact and vanishes if 013 = a\2- The 
coefficient of the contact diverges in the limit a,\s — > 0, 
which indicates that the mean shift in the rf frequency is 
sensitive to the range tq of the interactions between the 
atoms if (Z13 is not large. The sensitivity to the range 
arises in this case because r(w) has a high-frequency tail 
that decreases like w~ 3 / 2 and is proportional to the con- 
tact 0: 



11, 12] implies that the e + e annihilation into hadrons can be expressed formally 



r(«) 



n 2 



C\2 



(013 « 0) 



(2) 



The power-law tail extends out to oj of order l/(m?o), 
beyond which it is cut off by range effects. 

In this paper, we point out that the short-time oper- 
ator product expansion of quantum field theory provides 
a general method for deriving connections between few- 
body physics and many-body physics. We use it to red- 
erive and extend the universal results for rf spectroscopy 
in Eqs. (Ilb|) and ©. We also use it to derive new sum 
rules that are insensitive to range effects. 

The operator product expansion (OPE) is a power- 
ful tool for studying strongly-interacting quantum field 
theories that was invented independently by Wilson, 
Kadanoff, and Polyakov in 1969 14|. It expresses the 



product of local operators at different space-time points 
as an expansion in local operators with coefficients that 
are functions of the separation of the operators: 

O a (R +\r,T+ \t) O b (R — |r,T — \t) 

= Y,W c {r,t)O c {R,T). (3) 

c 

The sum is over the infinitely many local operators Oq- 
The functions Wc(r,t) are called Wilson coefficients. 
The short- distance OPE, with the operators at the same 
time (t = 0), is an asymptotic expansion for small \r\ 
15]. It can be generalized to a short-time OPE by ana- 
lytically continuing the time difference t to a Euclidean 
time r defined by t = —it. This form of the OPE is an 
asymptotic expansion for small \r\ and r. 

A classic application of the short-time OPE in high 
energy physics is electron-positron (e + e~) annihilation. 
In high-energy e + e~ collisions, most of the annihilation 
cross section is into hadrons (mesons and baryons). The 
fundamental theory for hadrons is quantum chromody- 
namics (QCD), which is a strongly- interacting quantum 
field theory. As far as the hadrons are concerned, the 
annihilation of an e + e~ pair with center-of-mass energy 
E results in the electromagnetic current operator 3(r,t) 
acting on the QCD vacuum state. This creates a quark 



in terms of the expectation value in the QCD vacuum of 
a product of electromagnetic currents: 



a(E) = (47ra/3£ 4 ) Imijdt e l{E+lt)t Jd 3 



x (0\T J (r,t) ■ J(O,0)|0), 



(4) 



where a m 1/137 is the fine structure constant of QED. 
The short-time OPE can be used to expand the opera- 
tor product J(r,t) ■ J(0,0) in powers of local operators 
with increasing dimensions. The leading operator is the 
unit operator, whose expectation value in any normalized 
state is 1 and whose dimension is 0. The next higher di- 
mension operators that have nonzero expectation values 
in the QCD vacuum are scalar quark operators with di- 
mension 3 and the gluon field strength operator with di- 
mension 4. Because of the asymptotic freedom of QCD, 
the Wilson coefficients can be calculated using perturba- 
tion theory in the QCD coupling constant a s . For large 
complex E, the Fourier transforms of the Wilson coeffi- 
cients of the higher dimension operators are suppressed 
by powers of E. The cross section at asymptotically large 
energy E can therefore be obtained by truncating the 
OPE after the unit operator and calculating its Wilson 
coefficient to leading order in a s . Inserting the OPE into 
Eq. (0} and neglecting quark masses relative to E, one 
obtains the simple result 



a(E) 



Ana' 



E 



(5) 



where the sum is over the quark flavors q and e q is the 
electric charge (+| or — |) of the quark [l6|. The leading 
corrections at large E decrease as powers of 1/ hi(E) and 
can be calculated using perturbation theory in a s . 

To apply the OPE to ultracold atoms, the problem 
must be formulated in terms of a local quantum field 
theory. Fermionic atoms with hyperfine spin states 1, 
2, and 3 that interact with large pair scattering lengths 
ai2, 023, and 013 can be described by a local quantum 
field theory with quantum fields ip cr (r) 1 a = 1,2,3. The 
interaction term in the Hamiltonian density that gives 
the large scattering length 012 is 

"Hint = (Ao,i2/m)V>Iv4 ? /'2?/'i- (6) 

The bare coupling constant Ao,i2 depends on the ultra- 
violet momentum cutoff A: Ao,i2 = 47r/(l/ai2 — 2A/V). 
Similar interaction terms give the large scattering lengths 
013 and a 2 3. In Ref. the contact Cyx associated with 
atoms 1 and 2 was identified as the integral over space of 
the expectation value of a local operator: 



Cn= d 3 R(\l 12 ^l4i; 2 MR))- 



(7) 
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Note that the contact density operator has an additional 
factor of the bare coupling constant compared to the in- 
teraction term in Eq. ([5]). In Ref. js|, two of Tan's uni- 
versal relations were rederived using the short-distance 
OPE. The OPE for tp\ and ip a was used to show that the 
momentum distributions for atoms of types 1 and 2 have 
identical l/k 4 tails whose coefficient is the contact C\2- 
The OPE for an( A Wi^>i was use d to show that the 

contact density is proportional to the number of pairs of 
atoms of types 1 and 2 per (volume) 4 / 3 . 

The universal relations for rf spectroscopy in Eqs. (jlbp 
and ([2]) can be rederived and extended by using the short- 
time OPE. We consider a many-body system consisting 
of ultracold atoms of types 1 and 2 only. As far as the 
atoms are concerned, rf spectroscopy with angular fre- 
quency ui proceeds through the action of the local oper- 
ator ipl^ir , t) on the many-body state. The operator 
transforms an atom of type 2 into an atom of type 3 
with energy larger by oj. The inclusive rate T(co) for pro- 
ducing final states containing an atom of type 3 can be 
expressed formally in terms of the expectation value in 
the many-body state of a product of operators: 

r(w) = n 2 Imifdt e l(u+ie)t [d 3 R [d 3 r 



x (T^ (R+±r, t) V>3>2 (R-hr, 0)) • (8) 

The Wilson coefficients in the OPE arc determined by 
few-body physics and can be calculated diagrammatically 
using methods described in Ref. 1 7j . The lowest dimen- 
sion operators that have nonzero expectation values in a 
state that contains no atoms of type 3 are 1/4^2, which 
has scaling dimension 3, and the contact density opera- 
tor Aq i2V'i'02' ( /'2'0i which has scaling dimension 4. The 
scaling dimensions of these and other operators can be 
deduced from the o pera tor-state correspondence of con- 
formal field theory [18] . The Fourier transforms of the 
corresponding terms in the OPE are 

Jdt e iujt Jd 3 r |r,t) ^ 2 (i?- |r,0) 

= (i/u) 4>lMR) + iWuiu) A^M^iC-R) 
+ .... (9) 

The Wilson coefficients can be determined by matching 
the matrix elements of both sides between asymptotic in- 
coming and outgoing few-atom states. For f/'IV^, one can 
use single-atom states of type 2. For the contact density 
operator, one can use two-atom scattering states consist- 
ing of atoms of types 1 and 2. The matrix elements for 
these states can be calculated using the analytic solution 
to the 2-body problem. The matrix element of the con- 
tact density operator is proportional to l/|l/ai2 + ik\ 2 , 
where k 2 /m is the total energy of the scattering atoms. 
Its Wilson coefficient is 



The OPE in Eq. © is an asymptotic expansion for large 
\u>\ along any ray in the complex u plane except possibly 
along the real u> axis, where the Wilson coefficients have 
poles and branch points. 

The scaling of the Wilson coefficients at large w is de- 
termined simply by the scaling dimensions of the opera- 
tors and by symmetries. The rf transition operators on 
the left side of the OPE in Eq. (|9|) each have scaling di- 
mension 3. Since the integral over nonrelativistic space 
and time has dimension —5, the total scaling dimension is 
1. On the right side of the OPE, the Wilson coefficient of 
an operator of dimension d must have terms that scale as 
UJ (i-d-n)/2^ a -iyn ^ w h ere a j s a scattering length. Unless 

factors of a -1 are required by a symmetry, the leading 
behavior is t i / 1 ~ d '/ 2 . For ip\ip2, which has scaling dimen- 
sion 3, the leading behavior is correctly predicted to be 
uj^ 1 . The contact density operator has scaling dimension 
4, so the leading behavior might be expected to be uj~ 3/>2 . 
However the rf transition operators are invariant under a 
subgroup of an SU (2) symmetry that is respected by the 
Hamiltonian if a\2 — ai3 10]. The contact density opera- 
tor does not respect this symmetry, so Wi2(w) must have 
a factor of a^ 1 — a^ 1 that vanishes when ai2 = 013. Thus 
its leading behavior is uj~ 2 , in agreement with Eq. (fT0|) . 

The short-time OPE can be used to derive sum 
rules for integrals over the rf frequency of the form 
J^^du! F(u))F(u>), where F(ui) is a weight function that is 
analytic on the real ui axis. Upon inserting the expression 
in Eq. (J8J) for T(ui) as a discontinuity in ui into the inte- 
gral, it can be expressed as a line integral over a contour 
that wraps around the real axis. The sum rules in Eq. ([T]) 
can be derived by deforming the contour into a circle at 
infinity in the complex oj plane. The short-time OPE 
converges everywhere along the contour except possibly 
near the real axis. Inserting the OPE in Eq. ([9]) into the 
expression for r(w) in Eq. ([8|) and then evaluating the 
contour integrals, we obtain the sum rules in Eq. (JT]). 

The short-time OPE can also be used to derive the 
large-frequency behavior of the rf intensity. Inserting the 
OPE in Eq. © into the expression for T(lj) in Eq. |8]), 
we obtain 



VL 2 {a 



12 



47ry / mw 3 / 2 (a 13 + muj) 



(11) 



'13 



\J — mijj 



Airmcu 2 a 13 - y/—muj 



(10) 



It was pointed out in Ref. [19| that T(cj) should decrease 
asymptotically like w -5 / 2 . The coefficient of the contact 
in the w~ 5 / 2 tail is a new result. The analytic result for 
the rf transition rate T(lj) for the weakly-bound dimer as- 
sociated with large positive ai2 agrees with Eq. (fTTj) 
in the limit to l/(ma 2 2 ) if we use the fact that the 
contact of the dimer is C12 = 87r/ai2 [7[. The result in 
Eq. (fTTj) for the tail in the rf transition rate relies on the 
scattering length 013 being large compared to the range 
ro. If ai3 is not large, the corresponding result can be 
obtained by taking the limit 013 — > in Eq. dill) , which 
gives Eq. ©. 
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Sum rules that are less sensitive to range effects can be 
obtained by using a weight function F(ui) that decreases 
as uj — > oo. Using a Lorentzian centered at ujq with half- 
width at half-maximum 7 gives a 2-parameter family of 
sum rules. Upon deforming the contour into a circle at 
infinity along which the integral vanishes, we are left with 
the contributions from the poles in uj at uo i *1- The 
convergence of the expansion for the sum rule is governed 
by the convergence of the OPE at this complex frequency 
and is therefore insensitive to the breakdown of the OPE 
near the real axis. In the limit a\3 — > 0, this sum rule is 



f 

j — { 



dui- 



7/tt 



(uj 



w ) 2 + 7 2 



PM = ^N: 



n 2 [{tu 2 - 7 2 )6+ + 2a; o7 b- - 2u 7a 12 1 ] 
Airmuj^ 



+ .... (12) 

where w 1 — (uj 2 + 7 2 ) 1 / 2 and b± = [m(u> 1 ± uj ) /2] 1 / 2 . 
The analogous sum rule fo_r_ e + e~ annihilation into 
hadrons was derived in Ref. 2l|. Sum rules of the form 
j£°duj F(uj)T(uj), where F( uj) is a polynomial, are also 
insensitive to range effects if uio <C l/(mrg). The con- 
vergence of the expansion for the sum rule is governed 
by the convergence of the OPE on the circle \uj\ = ujq. 
Sensitivity to the breakdown of the OPE near the real 
axis can be decreased by including a factor of uj — ujq in 
F(uj). 

In the OPE in Eq. ©, there are additional terms 
that have nonzero matrix elements in states that contain 
atoms of type 3. The lowest dimension operator is the 
23 analog of the 12 contact density operator in Eq. ((7|). 
If the rf operators in Eq. (|9|) are replaced by their hermi- 
tian conjugates, the local operators include V'lV'a an d the 
13 contact density operator. There is also a 123 contact 
density operator proportional to ^V'lV'jiV'sV^V'ij whose 
scaling behavior is governed by Efimov physics [ij. It 
has a complex scaling dimension whose real part is 5 
and whose imaginary part is ±2so, where sq « 1.00624. 
Thus its Wilson coefficient will decrease asymptotically 
like 



multiplied by a log-periodic function of u> with 
a discrete scaling factor of approximately 515. The cal- 
culation of its Wilson coefficient requires the numerical 
solution to a 3-body problem. 

The short-time OPE can be applied to other oper- 
ators to derive sum rules and high-frequency tails for 
the associated spectral functions. The OPE for t/jJ. and 
ip a provides universal information about the spectrum of 
single-particle excitations. The OPE for density opera- 
tors VctVv provides universal information about the spec- 
trum of density fluctuations. The OPE for two stress ten- 
sors provides constraints on the spectral functions that 
determine transport coefficients, such as the viscosity. 

Universal relations involving the contact provide non- 
trivial examples of aspects of many-body physics that are 
controlled by few-body physics. The OPE reveals these 



aspects by expressing observables in terms of Wilson co- 
efficients that are determined by few-body physics. Of 
the known universal relations, most have been derived 
from the analytic solution to the 2-body problem, but 
one has been derived from the solution to the 3-body 
problem j2^|. The derivation of new universal relations 
from numerical solutions to the 3-body and higher few- 
body problems presents an interesting challenge. 

Note added: The contact has recently been measured 
in experiments with ultracold 6 Li atoms [23j and 40 K 
atoms [24J. Several of the universal relations involving 
the contact have been verified experimentally. 

This research was supported in part by the Army Re- 
search Office and the Air Force Office for Scientific Re- 
search and by the Department of Energy. 



[2] 
[3] 
K 
[5 

[6] 

[r] 

[8 
[9 



[10 
[11 

[12 
[13 
[14 



[is; 
[ie; 
[17; 
[is; 

[19 
[20 
[21 



[22 
[23 
[24 



[1] E. Braaten and H. W. Hammer, Phys. Rept. 428, 259 

(2006) . 

C.A. Regal et al., Nature 424, 47 (2003). 
C. Chin et al., Science 305, 1128 (2004). 
Y. Shin et al., Phys. Rev. Lett. 99, 090403 (2007) . 
J.T. Stewart, J. P. Gaebler, and D.S. Jin, Nature 454, 
744 (2008). 

Q. Chen et al., Rep. Prog. Phys. 72, 122501 (2009). 
S. Tan, Annals of Physics 323, 2952 (2008); ibid. 323, 
2971 (2008); ibid. 323, 2987 (2008). 

E. Braaten and L. Platter, Phys. Rev. Lett. 100, 205301 
(2008). 

F. Werner, L. Tarruell, and Y. Castin, Eur. Phys. J. B 68, 
401 (2009); S. Zhang and A.J. Leggett, Phys. Rev. A 79, 
023601 (2009); R. Combescot, F. Alzetto, and X. Ley- 
ronas, Phys. Rev. A 79, 053640 (2009). 
Z. Yu and G. Baym, Phys. Rev. A 73, 063601 (2006). 
M. Punk and W. Zwerger, Phys. Rev. Lett. 99, 170404 

(2007) . 

G. Baym et al., Phys. Rev. Lett. 99, 190407 (2007). 
W. Schneider and M. Randeria, |arXiv:0910.2693] . 
K.G. Wilson, Phys. Rev. 179, 1499 (1969); 
L.P. Kadanoff, Phys. Rev. Lett. 23, 1430 (1969); 
A.M. Polyakov, Zh. Eksp. Teor. Fiz. 57, 271 (1969). 
K.G. Wilson and W. Zimmerman, Commun. Math. Phys. 
24, 87 (1972). 

T. Appelquist and H. Georgi, Phys. Rev. D 8, 4000 
(1973); A. Zee, Phys. Rev. D 8, 4038 (1973). 
E. Braaten, D. Kang and L. Platter, Phys. Rev. A 78, 
053606 (2008). 

Y. Nishida and D. T. Son, Phys. Rev. D 76, 086004 
(2007). 

A. Perali, P. Pieri, and G.C. Strinati, Nature Physics 5, 
736 (2009). 

C. Chin and P.S. Julienne, Phys. Rev. A 71, 012713 
(2005). 

E. C. Poggio, H. R. Quinn and S. Weinberg, Phys. Rev. 
D 13, 1958 (1976). * 

F. Werner a nd Y. Castin, |arXiv:1001.0774| . 

H. Hu et aL larXiv:1001.3200l 
J. T. Stewart et al, |arXiv:1002.1987 



